Experimental quantum state tomography of a solid state qubit 
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Full quantum state tomography is used to characterize the state of an ensemble based qubit imple- 
mented through two hyperfine levels in Pr 3+ ions, doped into a Y2SiOs crystal. We experimentally 
verify that single-qubit rotation errors due to inhomogeneities of the ensemble can be suppressed us- 
ing the Roos-M0lmer dark state scheme [lj. Fidelities above > 90%, presumably limited by excited 
state decoherence, were achieved. Although not explicitly taken care of in the Roos-M0lmer scheme, 
it appears that also decoherence due to inhomogeneous broadening on the hyperfine transition is 
largely suppressed. 

PACS numbers: 03.65.Wj, 03.67.Lx, 42.50.Md, 42.50.-p, 42.50.Dv 



A large variety of systems are presently investigated 
in order to find out whether they can be used as hard- 
ware for quantum computers. The present work is carried 
out on a solid state based system, rare earth ions doped 
into inorganic crystals. As in several other solid state 
systems the qubits are encoded in nuclear spin states, 
which for rare earths can have coherence times of seconds 
and where much longer coherence times are predicted [3] • 
For being a solid state system the rare earth ions are un- 
usual because their optical transitions can have coherence 
times as long as several ms |3j, |4| . Quantum state tomog- 
raphy have previously been carried out to characterize 
the fidelity by which superpositions on an optical tran- 
sition can be manipulated [5(. However, since coherence 
times for the hyperfine states are several orders of magni- 
tude longer, it is highly relevant to also investigate the fi- 
delity of arbitrary qubit rotations using hyperfine qubits. 
Multi-qubit gate operations can readily be implemented 
in the system, because optical excitation of an ion will 
induce frequency shifts >100 MHz (>10 4 line widths) of 
the optical transitions of nearby ions [f| . The large fre- 
quency shift of the optical transition makes it possible 
to entangle two nearby ions using operations with a du- 
ration of just a few ns Q- A scalable implementation 
of the rare earth ion scheme can e.g. be achieved using 
a short lifetime readout ion, acting as a state sensitive 
probe for the local environment [7| in a manner similar 
to how the electronic spin of an NV center can probe the 
nuclear spin states of surrounding C 13 ions However, 
because of the hour-long lifetimes of the rare earth spin 
states 0, QiJ , it is possible to also create qubits consisting 
of an ensemble of ions, all in a specific quantum state. 
Each such qubit can be selectively manipulated by opti- 



cal pulses [y, [12| . These ensemble qubits, which give 
strong readout signals, can be used to investigate general 
properties of the system. In this work ensemble qubits 



are used to experimentally carry out arbitrary rotations 
on the qubit Bloch sphere and the results are character- 
ized by full quantum state tomography. 
The relevant part of the Pr 3+ :Y2SiOs energy level dia- 
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FIG. 1: (color online) a) energy level diagram, b) qubit Bloch 
sphere, dark, \D), and bright states, \B), have been indicated, 
c) indicating the state transfer along paths 1 and 2 yielding 
an e* e phase shift of the bright state 



gram is shown in Fig. 1. The qubit states |0) or |1) are 
represented by two of the three ground state hyperfine 
levels and the qubits can also be optically excited to the 
|e) state, which has a lifetime of 164 us. Even if all ions 
in a specific qubit can be prepared in the |0) or |1) states 
using optical pumping, the ensemble approach brings ad- 
ditional problems because it has to be ascertained that 
all ions in a qubit have the same wave function. This 
is complicated by the facts that different ions within the 
qubit will; 1) have slightly different optical transition fre- 
quencies, leading to both different response to excitation 
pulses because some ions will be slightly off resonance 
and dephasing on the optical transition, 2) have slight 
differences in hyperfine transition frequencies, leading to 
different response to excitation pulses as well as dephas- 
ing on the qubit transition, 3) experience different optical 
field strength and have different Rabi frequencies, due to 
the spatial profile of the beam, which could result in dif- 
ferent ions experiencing different pulse areas, finally, 4) 
when implementing, e.g., two-qubit gates the interaction 
between an ion in the first qubit and the nearby ion in 
the second qubit will be different for different instances. 
However, the scheme for two-qubit gates is designed such 
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that it compensates for the inhomogencity of the ion-ion 
interaction [6j and complications 1 and 3 have been solved 
by employing the techniques described in Inhomo- 
geneities in optical transition frequency and differences 
in Rabi frequency are compensated for by using complex 
hyperbolic secant pulses (abbreviated as sechyp pulses 
through the rest of the text), which efficiently transfers 
states on the Bloch sphere from one pole to the other, 
provided the Rabi frequencies of the individual ions are 
above a certain lower limit. This was experimentally 
verified in Rcf [ll| and traces b) and c) in Fig. 2 also 



shows a transfer from the |0) state to the |1) state with 
about 96% transfer efficiency. However, to demonstrate 
arbitrary operations on the qubit Bloch sphere, which is 
the main objective here, it is also necessary to address 
complication 2), the inhomogeneous broadening on the 
qubit transition. Since the sechyp pulses only compen- 
sate for detunings in the resonance frequency and dif- 
ferences in Rabi frequency for pole-to-pole transfers on 
the Bloch sphere, Roos-M0lmer [l| introduced a basis 
change such that operations from an arbitrary point on 
the qubit Bloch sphere (Fig. lb) could be implemented 
as pole-to-pole transfers in the new base. The concept 
is schematically pictured in Fig. 1 and briefly described 
below. Two fields with optical frequencies ujq and oj\ and 
Rabi frequencies fig and f^i are driving the |0) — \e) and 
1) — |e) transitions, respectively. Adjusting the field am- 
plitudes such that £lo = fli creates one bright state, \B), 
and one dark state, \D) 

\B) = |0> - e-** |1) , \D) = |0) + |1) (1) 

depicted in the qubit Bloch sphere in Fig. lb. <f> is the 
relative phase difference between the two fields. The dark 
state wave function is not changed by the driving fields, 
however, the interaction between the driving field and the 
ions will drive the bright state along , e.g., path 1, on the 
\B) — |e) Bloch sphere (Fig. lc). If the fields Oo and f^i are 
sechyp pulses, they will compensate for detunings on the 
optical transition frequency for the \B) — > |e) transfer. If 
now the e) — ► \B) transfer is carried out along a different 
path on the Bloch sphere, path 2, separated an angle 9 
from path 1, the bright state will have undergone the 
operation \B) = e %e \B). In the qubit basis (|0), |1)), this 
is equivalent to the operation 



U 



cos I ie 1 ^ sin | 
sin % cos % 



(2) 



Thus we are able to carry out unitary rotations about 
any axis on the equator of the Bloch sphere representing 
the qubit basis. The angles 9 and <j> can be set explicitly 
in the experimental environment, giving us high accuracy 
in creating the qubit rotation. 

A Coherent 699-21 dye laser frequency stabilized against 
a spectral hole in a Pr 3+ :Y2SiOs crystal yielding a co- 
herence time >100 us and a frequency drift < 1 kHz/s 



was used for the experiments [13[ . The light passed twice 
through an 200 MHz acousto-optic modulator such that 
light pulses with arbitrarily chosen phase, amplitude and 
frequency patterns could be created without any spatial 
displacement of the output beam. These pulses were 
then sent through a 350 MHz AOM driven by two rf 
frequencies separated by 10.2 MHz, which is the split- 
ting between the qubit levels. In this way the first AOM 
produced the sechyp pulses and the second AOM dis- 
tributed this at frequencies and ui\, ascertaining that 
the pulses at the two frequencies had identical frequency 
chirp and amplitude variations (but different overall am- 
plitude, phase and center frequency). After the AOM the 
light was passed through a single-mode fiber, to clean up 
the spatial mode. A few percent of the light was split off 
after the fiber and used as a reference. The rest of the 
light (about 50 mW) was focused onto a 0.5 mm thick 
Y 2 Si05 crystal where 0.05% of the Y ions had been sub- 
stituted by Pr 3+ to a 1/e 2 spot diameter of ~100 u.m, 
yielding a Rabi frequency of maximum 2 MHz for the 
strongest transitions. The light transmitted through the 
crystal was imaged onto a 50 u.m pinhole only transmit- 
ting light from the center of the laser spot in the sample. 
Within this region the intensity varied by less than 30%. 
Reference and signal beams were detected by two Thor- 
labs PDB150A detectors and the signals from the two 
detectors were divided to reduce the effect of laser am- 
plitude fluctuations. The two-color sechyp pulses used in 
the experiment had a duration of 4.4 us and a FWHM 
of 1.2 us 



The qubits are created as in Refs [ill, [14|. Fig. 2a shows 
a part of the inhomogeneously broadened Pr 3+ absorp- 
tion line where all ions absorbing within an 18 MHz fre- 
quency interval have been transferred to other hypcrfme 
states through optical pumping by repeatedly scanning 
the laser back and forth in frequency. Ions within a nar- 
row frequency range are burnt back into the emptied fre- 
quency interval and placed in the |0) state, creating a 
peak with an inhomogeneous width of about 170 kHz 
(Fig. 2b). Because of the upper state hypcrfme splitting, 
the presence of ions in the |0) state shows up in the ab- 
sorption spectrum as three peaks separated by the upper 
state splittings, 4.6 MHz and 4.8 MHz. All spectra in 
Figs. 2 and 3 are recorded using the rapid chirp tech- 
niques developed in Refs [l5| and [lj| . Fig. 2c shows the 
absorption spectrum after first applying one sechyp pulse 
at the luq frequency to the |0) — |e) transition, bringing 
the entire population to the excited state, and then ap- 
plying a second sechyp pulse at the lo\ frequency on the 
|1) — |e) transition, transferring the population to the |1) 
state. The total transfer efficiency is about 96%. The 
homogeneous dephasing time, T%, for the optical transi- 
tions generally depends on the density of excited state 
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With a dopant concentration of 0.05 % the 
density of excited ions for one excited, 170 kHz wide 
qubit is about 3-10 14 /cm 3 . From Ref. [l8| this excited 



3 



1.5 
1.0 
0.5 



(a) 



-20 

0.6 
0.5 
0.4 
0.3 
0.2 

0.1 hf" 




(b) 



|l 2 



10.2 MHz_ 

17.3 MHz 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




1 0.2 MI[/_ 

17.3 MHz 



4 6 8 10 12 
Relative Frequency (MHz) 



FIG. 2: (color online) a) Tailoring of the absorption profile, 
all ions within a 18 MHz spectral interval have been removed 
by optical pumping, b) a qubit is created by returning ions, 
which have their |0) — |e) transition within 100 kHz from rela- 
tive frequency 0, to the optically pumped region, c) The qubit 
has been transferred from the |0) to the |1) state, the transfer 
efficiency is about 96% 
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FIG. 3: (color online) State tomography of the |0) — |1) state. 
Traces a), b) and c) show the projection of the state on the 
x-, y, and z-axes, respectively. Numbers 1, 2, 4 and 5 in trace 
c), refer to the numbering of the transitions in Fig. 2. Further 
explanation is given in the text. 



state density would give a excited state dephasing time 
of about 50 us, which is consistent with our own photon 
echo measurements of the dephasing time when eciting 
one qubit. Simulating the state-to-state transfers with 
the sechyp pulses used and a T2 of 50 us using a Bloch 
equation model, gives a maximum transfer efficiency of 
about 96%. It is consequently reasonable to assume that 
excited state dephasing is the main limiting factor for the 
transfer efficiency in Fig. 2c. 

We are now ready to use the two-color sechyp pulses, em- 
ploying the scheme outlined in Fig. 1 for performing arbi- 
trary single qubit operations on qubits like the one shown 
in Fig. 2b. Starting from state |0) (Fig. 2b) five different 
states on the Bloch sphere were prepared (|0) + |1)), 
(|0) - |1)), (|0) + i |1)) and (|0) - i |1)) (normalization 
factors have been omitted for notational simplicity). To 
characterize the created state we calculate the fidelity, 

F, as F = (*thoory| Poxp l^thoory), where ^theory is tllC 

desired state and p exp is given by [3| 

Pcxp = 0.5 • [tx{p)I+tv(Xp)X+tv{Yp)Y+tx{Zp)Z]. (3) 

I is the identity matrix, X, Y and Z are the Pauli matri- 
ces and tr(Xp), tr(Yp) and tr(Zp)) are the experimental 
results from measurements of the projection of the pre- 
pared qubit state on the x-, y- and z-axes and tr(p) is set 
equal to unity (It is here assumed that decay to the |aux) 
state, Fig. la, can be neglected and indeed the branching 
ratio to this state is small I n Fig- 3 the arrow in 

the direction of the negative x-axis in the Bloch sphere 
schematically illustrates that the (|0) — 11)) state has been 



prepared. Traces a), b) and c) shows the results from the 
measurements of the projections of the (|0) — |1)) state 
on the x-, y- and z-axes. The only experimental mea- 
surement at hand is to measure the projection of a state 
on the z-axis, i.e., a frequency resolved absorption mea- 
surement telling the fraction of ions in the |0) and the 
|1) states. To measure the state projection on the x- and 
y-axes it is then necessary to carry out rotations on the 
Bloch sphere projecting these axes onto the z-axis. How- 
ever, projecting the x-axis on the z-axis requires exactly 
the same pulses as creating the (|0) — |1)) state from the 
|0) state and projecting the y-axis on the z-axis requires 
exactly the same pulses as creating the (|0) — i |1}) state 
from the |0) state. Thus, there is now sufficient infor- 
mation to carry out the single-qubit operation and the 
state tomography. First a qubit is prepared as in Fig. 2a, 
then an operation on the \B) — \e) Bloch sphere is done 
following the procedure in Fig. lc where the angles 4> and 
9 are chosen according to Eq. [3] in order to carry out the 
desired operation. The qubit preparation and the qubit 
rotation is carried out three times. The first time the 
rotation is followed by a second rotation projecting the 
x-axis on the z-axis (followed by an absorption measure- 
ment, compare Fig. 3a), the second time the rotation is 
followed by a rotation projecting the y-axis on the z-axis, 
(compare Fig. 3b) and finally the rotation is followed by 
a projection of the z-axis on the z-axis by choosing 9 in 
Fig. lc equal to zero, bringing the bright state up and 
down along the same path. The z-axis projection can of 
course be measured by an absorption measurement di- 
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rectly, without any operations on the bright state, but 
actually carrying out the operation to project the z-axis 
onto the z-axis using = 0, means that the X, Y and Z 
tomography measurements are all carried out in an equiv- 
alent manner. However, using only the pulses above the 
fidelity would be very low. The action of the two-color 
sechyp pulses are to first promote the bright state to the 
excited state. In the general case the ions will then be in 
a superposition of the dark state and the excited state. 
Because of the inhomogeneous broadening of the opti- 
cal transition, different qubit ions will acquire different 
phase factors while being in this superposition. However, 
a phase factor can be considered as global, and thus dis- 
regarded, if it appears in front of all qubit states. After 
the excited state has been returned to the bright state 
along path 2 in Fig. lc , the detuning dependent phase 
factor accumulated in the excited state can now also be 
accumulated on the dark state. By sending in two new 
two-color sechyp pulses, identical with the previous pair 
except that, first, the phase <fi (see Eq.[T]) now is increased 
by an amount tt (this means the new pulses will act on the 
dark state), second, 9 is set equal to zero, which means 
the dark state is taken up and down on the \D) — \e) Bloch 
sphere along the same route. Consequently U in Eq. [2] 
will (apart from an overall phase factor) for this second 
sechyp pulse pair, be an identity operation. However, 
the dark state part of the wave function will still have 
acquired the excited state phase factor due to inhomo- 
geneous broadening. Dephasing due to inhomogeneous 
broadening on the optical transition is now eliminated. 
A further discussion of this issue can be found in []]]. 
Thus, each single-qubit rotation as well as each tomog- 
raphy operation consists of four two-color sechyp pulses. 
Two on the bright state and two on the dark state. In 
the present experiment each trace in Fig. 3a-c is the re- 
sult after applying eight two-color sechyp pulses, where 
each such pulse is 4.4 us long, the total sequence for each 
trace then is 35.2 us. The fidelities obtained are given in 
Table I. Two different fidelities are given Fqr+qst is the 



Fqr 
Fqr+qst 



|1) ]0) + |l) |o> — 1 1> |0)+*11) |0)-»|1) 



0.96(2) 0.93(1) 0.93(1) 0.92(2) 0.91(2) 
1.02(2) 0.92(3) 0.87(1) 0.87(2) 0.85(4) 0.84(4) 



TABLE I: Fidelities for the single qubit rotation, Fqr, and 
single qubit rotation + quantum state tomography, Fqr+qst 
for five different states. A quantum state tomography for the 
starting state, |0), gives unity fidelity (within one standard 
deviation). 

fidelity calculated according to Eq. O where indices, QR 
and QST stand for Qubit Rotation and Quantum State 
Tomography, respectively. However, since the QR+QST 
operation just is two consecutive QR operations, it is rea- 
sonable to state the fidelity for a single qubit rotation, 
Fqr, as (Fqr+qst) 1 / 2 - Thus the fidelities for the sin- 
gle qubit operations are estimated to lie between 0.9 and 



0.96, giving an average single operation fidelity of 0.93 
The fidelities obtained are remarkably good considering 
the dephasing times of the system. The average time 
spent in the upper state during the QR+QST opera- 
tion is t u = 8.8 us. e~*"/ T2 = 0.84, which, assuming 
p = 0.84p(prepared state) + 0.16p(mixed state), would 
give a best case fidelity of 0.92. Optically detected Free 
Induction Decay (FID) at the qubit transition, was used 
to determine the dephasing time on the hyperfine (hf) 
qubit transition, T2(hf). Using a two-color sechyp pulse 
the qubit was put in a superposition state. The coher- 
ence on the qubit transition was probed by a delayed 
optical probe pulse. The 10.2 MHz beat signal due to 
the coupling between the qubit coherence and the trans- 
mitted optical pulse was detected. The strength of the 
beat signal as function of probe pulse delay measures the 
dephasing on qubit transition. The complete QR+QST 
sequence is 35.2 us. From the FID measurement only 
about 20% of the qubit coherence remained after 35 us. 
The decoherence is reversible and caused by the inho- 
mogeneous broadening of the qubit transition. In view 
of the fact that only about 20% of the the qubit coher- 
ence remains after 35 us, the fidelities obtained may seem 
remarkable. Possibly, this could be explained by the re- 
sults in Ref, [2fj| , where it has been shown that dynamical 
Stark shift occurring during qubit rotation can suppress 
errors due to inhomogeneous shifts of the qubit levels by 
as much as a factor of 10. It could also be a Zeno ef- 
fect, which has been encountered previously in similar 
systems, such as NV centers 21(. Still it would definitely 



be interesting to further investigate the effect of the in- 
homogeneous broadening on the QR fidelity. 
The single qubit rotation fidelities could be improved; 1) 
by using pulses of shorter duration developed by optimal 
control theory [i^, [3], 2) by using the single instance 
scheme 0] which eliminates errors due to inhomogeneous 
broadening on the hyperfine and optical transitions, or 
3) by instead of the Pr ion use the Eu ion, where the up- 
per state dephasing time is an order of magnitude longer. 
Assuming the fidelities in this work are limited by T2 and 
T2(hf), one or several of these changes should enable fi- 
delities above 0.99. To get fidelities significantly beyond 
this value, harder focusing of the light beam to increase 
the Rabi frequencies and /or dark state schemes not pop- 
ulating the excited state 24 , 25| , would probably need to 
be used. 
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